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Abstract. We propose an alternative approach to the dissipative vortex dynamics 
occurring in a superfluid vortex lattice at finite temperatures. Focusing upon the 
pseudomomentum of a vortex and its surrounding quasiparticles, wc derive an equation 
of motion which, in spite of yielding the same evolution as the usual one for massless 
vortices, does not involve the vortex mass. This picture could provide further insights 
into the controversy about the nature of the vortex mass. 
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1. Introduction 

Below the lambda transition, at temperatures below T\ = 2.172 K, liquid helium is 
known as helium //and can be regarded as a mixture of a normal fluid with mass density 
p n and a superfluid with mass density p s . Such densities are temperature dependent, so 
that p s (T\) = Pn(0) = 0, while the total density p = p s + p n remains nearly constant. 
The most striking property of such a superfluid component is, perhaps, that it can only 
rotate through vortices having microscopic cores and quantized circulations. In fact, 
the circulation of the superfluid velocity field around each of such vortices is quantized 
in units of h/m^, the so-called quantum of circulation k, given by the ratio of Planck's 
constant and the mass of one 4 He atom. In practice, however, we only find configurations 
of one quantum per vortex, since they are favored by energetic considerations pQ. When 
a rotating sample of liquid helium is cooled below the lambda temperature, all the 
rotation of the superfluid becomes concentrated in such vortices, which eventually form 
a lattice consisting in a uniform array of vortex filaments parallel to the axis of rotation 
[H [2], [3] . However, the macroscopic superfluid velocity field, which corresponds to spatial 
averages over regions large compared with the spacing between vortices, yields the usual 
configuration of solid body flow, v s (r) = f2 rot z x r for a rotation frequency f2 rot around 
the z axis. On the other hand, the microscopic superfluid velocity field, ie without 
averaging, is irrotational except where a vortex is located. Since the circulation of such 
a field around a vortex yields the quantum of circulation k, according to the Stokes' 
theorem we may write 



where the area A in the x-y plane contains Na vortices and w denotes the microscopic 
vorticity, ie the curl of the microscopic superfluid velocity. On the left-hand side of (pQ) 
we have an averaged vorticity which should be identified with the macroscopic value 
V x v s = 20 rot z. Thus, the above equation yields a link between the macroscopic 
and the microscopic views of the superfluid component, namely the number of vortices 
per unit area corresponds to the ratio of the macroscopic vorticity to the quantum of 
circulation, Na/A = 20 rot //« (Feynman's rule [HE]). 

Just as the superfluid flow is microscopically formed by vortices, the normal fluid 
consists of superfluid quasiparticle excitations, phonons and rotons, the average flow of 
which is characterized by a normal fluid velocity field v n given by f2 rot z x r. That 
is, both fluids are expected to move with the same velocity at equilibrium. It is 
important to recall, however, that there are remarkable exceptions to this result, eg 
the equilibrium configuration of the Hess-Fairbank experiment [5], or the metastable 
superflow states having effectively infinite lifetimes pjj. Now, keeping the focus on 
the simpler situation of having identical normal and superfluid velocity fields at an 
equilibrium state to be reached within a finite relaxation time, one may interpret that 
such a behaviour arises, from a microscopic viewpoint, from the vortex motion with 
the normal fluid velocity in order to avoid dissipation. That is, any relative motion 
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of vortices with respect to the normal fluid in their vicinity, should be subjected to a 
friction force that causes such a motion to eventually cease. Such a mutual friction 
force [7j between the two fluids appears then as playing a central role in the mechanism 
which maintains the stability of the above equilibrium state. A phenomenological model 
for this macroscopic dynamics was proposed long ago through the so-called Hall-Vinen- 
Bekharevich-Khalatnikov (HVBK) equations [7J [8] . These basically consist of a Navier- 
Stokes equation for the normal fluid and an Euler equation for the superfluid, which, 
in the absence of pressure and temperature gradients, are coupled together only by 
a mutual friction term [H [2]. The original proposal of the HVBK model was later 
rederived from first principles within the framework of classical continuum mechanics 
[9], but a derivation from a full microscopic theory is still lacking. On the other hand, 
with the exception of a few works [TUl [TT] and because of their complexity, the HVBK 
equations have been mainly utilized so far to model helium II with a spatially uniform 
configuration of vortices. 

The simplest departure from the rotating equilibrium configuration is given by a 
normal fluid field of the form (f2 rot + w ) z x r, but unfortunately a simple ansatz like 
v n (r, t) = Q n (t)zxr does not constitute an acceptable solution of the HVBK equations. 
Simplicity, however may be preserved by changing to a different geometry consisting in 
rectilinear flows of uniform vorticity, 



(y < 0), where Q s (t) and fl n (t) should converge for t — > oo to a common steady 
state value. Then, to make contact with the standard rotational configuration, we may 
identify such a value with the former angular velocity f2 rot . Note that this assignment 
leads to a uniform vorticity V x v s = 2 Q Iot z, which coincides with that of the rotational 
scheme. 

The above macroscopic view of the interaction between superfluid and normal 
fluid, ruled by the HVBK equations, is complemented by the microscopic picture of 
a dissipative vortex dynamics arising from the scattering of quasiparticles by vortex 
lineal In the usual approach, the vortex equation of motion arises simply from assuming 
a vanishing total force, which is given by the sum of a hydrodynamic Magnus force and 
a dissipative force. Such a neglect of the vortex inertia corresponds to the assumption 
that its mass is given by the hydrodynamic mass of a core of atomic dimensions [I] . This 
result, however, has never been experimentally confirmed owing to the difficulties that 
embodies a direct measurement of the vortex mass [H]. Moreover, there are different 
theories [151 EE Ej that yield several orders of magnitude higher values for the vortex 



mass, casting doubt on models based on massless vortices^. On the other hand, it 
has been recently suggested that an unambiguous vortex mass may not exist, and that 

| Vortex bending in the form of thermal excitation of vortex oscillation modes, or collective excitations 

as Tkachenko waves are not expected to be relevant to this discussion [H [TJl [13] . 

§ There have also been conflicting results for the vortex mass in superconductors [18] . 
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inertial effects in vortex dynamics may be scenario- dependent [T5] . Given such an open 
debate, it seems to be quite advisable to follow an eclectic procedure, assuming in what 
follows a finite vortex mass per unit length, which we shall denote through the parameter 

77V 

There is a close analogy between the dynamics of massive rectilinear vortices and the 
well-known electro dynamical problem of point charges subjected to a uniform magnetic 
field and a perpendicular electric field. More precisely, there exists a whole mapping 
by which a 2-D homogeneous superfluid can be mapped onto a (2+l)-D electrodynamic 
system, with vortices and phonons playing the role of charges and photons, respectively 
[T5l 1201 121] 1221 123] . Here we shall only make a restricted use of this mapping, which 
corresponds to the formal analogy between the Magnus and Lorentz forces. Thus 
we may assume that the dissipative dynamics of the vortex lattice should be ruled 
by three characteristic frequencies, namely the imposed rotational frequency f2 rot , an 
initial departure Wq from this frequency and the cyclotron frequency stemming from the 
electromagnetic analogy [13] 

Vt = p s K/m v . (4) 

Taking as a lowest estimate of the vortex mass the value arising from the above 
hydrodynamical model, we have fl < 3 x 10 12 s _1 [lj. On the other hand, typical 
experimental values are of order Q rot ~ 1 s _1 , so we shall restrict our study to cyclotron 
frequency values fulfilling Q rot <C Q. 

The above assumption of massive vortices and rectilinear flows leads us to a very 
useful magnitude, which to the best of our knowledge has not been utilized so far, 
that is the concept of vortex pseudomomentum. In fact, such a pseudomomentum 
corresponds to the vortex generator of translation, as can be straightforwardly shown 
from the electromagnetic analogy [24J. More generally, the translation generated by a 
pseudomomentum corresponds to a motion of the physical state (vortex) but keeping 
the medium, the uniform superfluid in this case, fixed [25]. In addition, contrary to 
the other two momenta (canonical and dynamical) that can be ascribed to a vortex, 
the vortex pseudomomentum turns out to be free from the ambiguities carried by the 
vortex mass. Here it is worthwhile noticing that a similar situation occurs in the case of 
the normal fluid. In fact, it can be shown that the momentum of a sound wave (or of a 
phonon in quantum mechanics), turns out to be a very complicated object, which may 
not even have a well-defined value at all, while its pseudomomentum is a simple quantity 
and far more useful [25], [26] [27] . This led us to investigate a pseudomomentum approach 
to the dissipative vortex dynamics occurring in a uniform vortex lattice, finding that 
this formalism leads to the same evolution as that predicted by the usual approach, 
whereas it involves far less restrictive assumptions about the value of the vortex mass. 

This paper is organized as follows. In the next section we study a solution of 
the HVBK equations for the above rectilinear flows. In section [3] we analyze the 
dissipative vortex dynamics from the viewpoint of the usual phenomenological approach. 
A Hamiltonian approach is proposed in section H] which consists of a vortex Hamiltonian 
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based on the analogy to electrodynamics (section 14. ip and a Hamiltonian for the 
quasiparticle gas representing the normal fluid (section [4.2l) . Finally a pseudomomentum 
equation of motion is obtained and discussed in section 14.31 



2. Two-fluid equations 

The HVBK equations for the rectilinear flows ([2])- ([3]) are very simple and read 

Pn^ = F (5) 

Ps~^ = ~ F, (6) 
where the mutual friction force [H [21 [7], 

p _ PnPs_ 0sW[ _ B(Vn - Vs ) + B 'zx (v„ - v s )], (7) 

becomes the single coupling term between both fluids and the dimensionless coefficients 
B and B' weight dissipative and nondissipative contributions to such a force. If we 
assume that the system is not far from equilibrium, so that the factor Q s (t) in ([7j) can 
be approximated by O ro t , replacing in (JSJ) to (JTj) the velocity fields according to (j2D and 
([3]), an elementary calculation leads to the solution 

a n (t) = n vot + w e~ n ^ Bt (8) 

Q s (t) =Q TOt -^w e- n ^ Bt , (9) 
P 

where we should assume temperatures below 1 K, so that the normal fluid density turns 
out to be much smaller than the superfluid one (p n / p s < 1CT 3 ) and the last term in (J9j) 
becomes negligible, ie 

^ « I- (10) 

Note that this model of strictly rectilinear flows amounts to ignoring nondissipative 
contributions to the mutual friction force, ie it leads to B' = 0. 



3. Dissipative vortex dynamics: phenomenological approach 

A microscopic view to the above lattice leads us to focusing on a single vortex dynamics. 
Thus, we assume that the motion of each vortex is ruled by p]: 

m v r = p s KZ x (f - v a ) - D(y - v„), (11) 

where r = ix + yy denotes the position of the vortex core and D denotes a microscopic 
friction coefficient arising from vortex-quasiparticle scattering. The right-hand side of 
(TTTj) represents the total force acting on the vortex, namely the Magnus force (first 
term) plus the dissipative force (second term). Here we have disregarded again any 
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contribution from a nondissipative component of the mutual friction forcqjJJ. Now, 
replacing ©, © and fi s (t) ~ fi rot in ([Til we obtain, 

x = - fi— x -Sly -2 a-ot ^— y ( 1 + ^- e - Q ^ Bt ] (12) 

y = -ft— y + 2 fi rot ?/ + x. (13) 
The equation (fT2l) can be formally solved in x(t), 
x(t) = e- {nD/p ° K)t x(0) - f dre- inD/p ° K)T 



x 



Q y(t - t) + 2 a o t n— y(t -r)(l + ^e"^ 8 ^)) 



(14) 



and replacing the above expression in ([13]) we are led to a second-order integro- 
differential equation in y(t). This problem, however, may be greatly simplified by 
utilizing the Markov approximation. In fact, we first note that the coefficients in the 
arguments of the exponentials in (fT4"|) define two very different time scales, namely a 
microscopic one given by (QD / p s n)~ l and a macroscopic one given by (VL mt B)~ l , ie 

« 1, (15) 

where we have taken into account that the macroscopic and microscopic friction 
coefficients are related for T< 1 K through [T], 

B = —. (16) 

Note that the inequality ffT5l) is expected to break down at extremely low temperatures. 
Now, assuming that t in (|14p belongs to the macroscopic time scale, we have that 
the exponential factor in the first term will be negligible, while the same exponential 
function, which depends on r in front of the integrand, will make that only the lowest 
portion of the integration domain (r <C t) gives a nonnegligible contribution to the 
integral. That is, we may safely approximate all the dependencies on t — r in (|14|) by t, 
and the upper limit of the integral by +oo. Finally, the Markov approximation to such 
an equation reads as, 



x(t) = - / dTe-^ Dtp ° K)l 



Qy(t)+2Q TOt Q—y(t) (i + ^ e -^Bt 

p s K \ il rot 



y-2n vot y l + _JL e -»«*"t ; (i 7 ) 



D V f^rot 

which replaced in f|T3|) and using (JTBj) yields, 

^y + y+*2Jt« be -<W« y = 0. (18) 
2p s il p s 

\\ Actually, according to some theories [28] and related experimental evidence [29l [30], such a 
component may in fact be non-existent. 



Vortex pseudomomentum and dissipation 



7 



Here it is useful to change to the adimensional macroscopic time variable T = Q ro tB t, 

and recall that p n /p s ^ 10 3 and B < 1 for T < 1 K. Then, assuming u> < ^ ro t and 
taking into account (Tl0|) and (1T51) . we realize that we may safely drop the term containing 
the second derivative in (1191) . since we may estimate that the coefficient in front of such 
a derivative will be a quantity of third order in comparison to the first-order small 
parameter in front of the exponential. Note that according to (jlj), this approximation 
turns out to be equivalent to the usual one of neglecting the vortex mass. Thus, we are 
led to a first-order differential equation which is easily integrated yielding, 

y(t) = y(0) exp[( Pn w / Ps n iot )(e~ n ^ Bt - 1)], (20) 

ie a small vortex displacement in the direction perpendicular to the velocity of the 
background superflow, while ( iTTjl becomes 

x = -2Q Iot y, (21) 

which simply states that the x-component of the vortex velocity will coincide with the 
superfluid velocity. 



4. Dissipative vortex dynamics: Hamiltonian approach 

4-1. Vortex Hamiltonian 

We start from the vortex equation of motion (fTTl) in the absence of normal fluid 

r = fizx(r-v s ), (22) 

which turns out to be analogous to that ruling the two-dimensional motion of a negative 
point charge in the presence of magnetic and electric fields in the z and y directions, 
respectively [23]. Such an equation derives from the Hamiltonian 

H v = -y (i£ + v 2 y ) - n iotPs Ky 2 , (23) 

being (Landau gauge) 

Px n 
v x = lly 

Vy = -S 

m v 

where p = p x ~k + p y y corresponds to the vortex canonical momentum. Note that such a 
momentum as well as the Hamitonian (1231) are given per unit length of the vortex line. 
Then, from Hamilton equations it is easy to check that the expressions (1241) correspond 
to the velocity of the vortex core r, while the acceleration is indeed given by ( 122]) . From 
the above electromagnetic analogy it is also useful to represent the coordinate r as the 
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sum of the center coordinate r = £o x + y Q y of the cyclotron circle plus the relative 
coordinate r' from such a center [21], being 

Vy 

v x Q (25) 

Then, it is easy to extract the time evolution of such coordinates working to zero-th 
order in the small parameter Q TOt /Q: 

r = -2fi rot ?/ x, (26) 

which means that the center coordinate r will move with the superffuid velocity, while 
r'(t) will perform a counterclockwise circular cyclotron motion with angular frequency 
Q. Finally, we may see from (125]) that the limit of a vanishing vortex mass (Q — ► oo), 
which is often found in the literature [31], amounts to ignoring the cyclotron motion 
(r' -> 0, r = r ). 



4-2. Quasiparticle Hamiltonian 

Although we have described the vortex Hamiltonian in classical terms, it is immediate 
by means of the electromagnetic analogy to switch to a quantum mechanical picture 
[21]. This constitutes a necessary generalization, since the normal fluid will be treated 
as a low-temperature quantum gas. In fact, such a fluid will be represented by the 
following Hamiltonian: 

H n = ^ a q a q' ( 27 ) 
q 

where a q denotes a creation operator of quasiparticle excitations of (pseudo-) momentum 
hq and frequency ujq. Measuring such a frequency from the lab frame, it is written as a 
Doppler-shifted frequency from the superffuid frame, Uq = u q + q • v s , where u q is the 
familiar (isotropic) dispersion relationship of 4 He quasiparticle excitations. Note also 
that we disregard any interaction between the quasiparticles themselves, since we shall 
work at low enough temperature, so that they remain dilute allowing their treatment 
as a noninteracting gas. 



4-3. Equation of motion in terms of pseudomomentum 

Note that according to the electromagnetic analogy, there are three kinds of momentum 
to be ascribed to the vortex, viz the canonical one p, the dynamical one m v v = p s nzx r' 
and the so-called pseudomomentum [21], which is given by K = —p s nz x r and should 
be regarded as the generator of translation. Then, adding such a pseudomomentum to 
the quasiparticle pseudomomentum ^2 a q a q; we have the pseudomomentum of the 
whole system. Note that only the x-component of the vortex pseudomomentum p s Kyo 
will commute with the vortex Hamiltonian ( 1231) . unless Q rot = 0. This result may be 
easily interpreted, since a superflow of velocity v s = —2 Q TOt y x produces a translation 
symmetry breaking in the ^/-direction. 
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Next let us analize the conservation of the x-component of the pseudomomentum 
for a given vortex of the lattice and its surrounding quasiparticles^J 

p s KL(y ) + n q h(\ • x = const., (28) 
q 

where n q = (at a q ) denotes the average number of quasiparticles with pseudomomentum 
hq and L denotes the vortex line length. We shall assume that such a population is well 
described by a local equilibrium form: 



4fc B Tsinh 2 (^; ? /2fc B T) 

where the quasiparticle frequency is measured from a reference frame where the local 
normal fluid is at rest, cu q = oj q — q- [v n — v s ], with v s = —2 f2 rot y x and v n = —2 f2 n (t) y x, 
[fl n (t) given by ©]. Then, replacing ([29]) in ([28]) we have, 

p s n(yo) + A v p„(v„, - v s ) • x = const., (30) 

where 

Pn = 12A v Lk B T S sinh^^^^T) ^ 31 ^ 

yields the normal fluid density and A„ = (2fi rot /fi;)~ 1 corresponds to the area per vortex 
of the lattice (Feynman's rule). Taking into account that the velocity fields in ( l30i) must 
be evaluated at y — (y ), and using the solutions (jHJ) and © we obtain, 

(Vo(t)) = (l/ (0)> [1 + (p„^o/p s ^ot)(e-^ tBt - 1)], (32) 

which according to (1101) . turns out to be equivalent to (1201) . In other words, starting from 
the conservation of pseudomomentum for a vortex and its surrounding quasiparticles, 
without any further assumptions regarding the vortex mass, we have shown that the 
y-component of r presents the same dissipative evolution as that obtained in section [3] 
for a vortex of negligible mass, for which r = r . Thus, we may see how the dynamics 
in terms of the pseudomomentum makes the coordinate ro the focus of attention, 
while the less relevant cyclotron dynamics of r', which was forced to vanish in the 
phenomenological treatment through the neglect of the vortex mass, is now naturally 
ignored. 

To study the dynamics along the x-direction we must focus on the y-component 
of pseudomomentum, which is given by the nonconserved vortex contribution —p s n(xo) 
(cf (1261) ). Finally, taking into account the whole vector pseudomomentum K tot of the 
vortex plus surrounding quasiparticles, we may write the following equation: 

= Fext ' (33) 

Here it is worth comparing with a recently derived conservation theorem for wave pseudomomentum 
plus a suitably defined classical vortex impulse [27] , 
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where 

he 

K to t = -p s Kz x (r ) + -tt— p n (v n - v s ) (34) 

'rot 

and 

F ex t = -Ps^z X v s (35) 

represents the force exerted on the vortex by the superfluid current [3], which is 
responsible for the nonconserved ^/-component of pseudomomentum. 

The equation of motion fl33|) constitutes the central result of this paper, which 
may be regarded as arising from a straightforward combination of the electromagnetic 
analogy and the fruitful concept of pseudomomentum. Here it is important to remark 
that contrary to the phenomenological equation of motion (TTTT) . the pseudomomentum 
equation (1331) does not involve the vortex mass, a fact to be regarded as most welcome, 
given its apparent elusive and ambiguous nature [19J. 

To conclude we may observe that it would be important to generalize this picture 
for the rotating superfluid, where the key magnitude would be represented by an angular 
pseudomomentum. To achieve this goal, given the drawbacks of the HVBK equations 
pointed out in section [TJ it would be helpful previously to seek for a validation of the 
present results without utilizing such equations as a starting point. 
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